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Abstract 

The model system manifesting phenomena peculiar to complex analytic maps is of- 
fered. The system is a non-autonomous ring cavity with nonlinear elements and filters, 
Keywords:M&ndelbrot and Julia sets; complex analytic maps; ring systems. 

When one studies simple quadratic map of a complex variable and with complex parameter 

z n+ i = c + z 2 n . (1) 

the interesting phenomena can be observed. The complex phase plane of a variable z, appears 
subdivided to two domains. The start from one area gives trajectory z n , escaping to infinity. 
The start from another area gives trajectory, remaining restricted: the trajectory wanders not 
far away from an origin. This second domain is named a Julia set. In a trivial case with c = 
it looks as a unit disk. With nonzero parameter values the Julia set is, fractal with non-smouth 
border, and has self-similar structure (see Fig. 1). 

The restricted in a phase space dynamics is possible not with any parameter c values. More- 
over, with different parameter values it has different properties - it can be periodic or chaotic. 
The useful method for examination of a system is the drowing of a chart of a complex parameter 
plane, for example as a chart of dynamical regimes. At Fig. 1 such chart is represented. On 
the plane c the object known as a Mandelbrot set is arises. It represents a fractal domain, 
parameter values from which corresponds to restricted dynamics. Mandelbrot set looks as a 
cardioid with a set of the attached round lobes (corresponded to periodic dynamics). This 
"cartus" is enclosed by fractal "mane", corresponded to chaotic dynamics. 

The huge amount of the mathematical literature is devoted to Mandelbrot and Julia sets 
[1-3]. The important task for physicist is the development of physical applications of the theory 
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Figure 1: Mandelbrot set (left panel) and Julia sets (right panels) for the complex quadratic 
map (1) at parameter and phase plane charts. Yellow areas correspond to the periodic dynamics 
of the complex variable z in a bounded domain (the periods are marked by figures). Black color 
designates bounded chaotic dynamics. Gray colors correspond to the escape of the trajectories 
to infinity by different steps of time (more dark color corresponds to more long escape time) 
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Figure 2: Diagram of the ring device with two nonlinear elements (violet rectangles of length Ri 
and R2), two filters, absorbing the signal of frequency 2u (blue box) and the signal of frequency 
oj (pink box). C and D are the external pumping signals of frequency u. Z n is a complex slow 
amplitude of the signal (of frequrncy u>) near after the knot, marked as 1 

of complex analytic dynamics and building up of the physical devices and systems in which the 
complex dynamics is implemented. Just several examples of such systems are exists [4-9] . In 
the present work a new example of a system in which Mandelbrot set can be implemented is 
offered. It is a non-autonomous ring system with nonlinear elements and filters. This model 
can be constructed for example on the basis of an optical system of the Ikeda cavity type [10] 
or as a ring system with the ferromagnetic structures as nonlinear elements [11]. 

Let us consider the system represented at Fig. 2. The signal can travel between knots 1- 
2-3-4. In the optical resonator these are the mirrors, turned around from each other so that a 
ray of light transits on a ring path. An external pumping C and D exciting the signal in the 
system are brought at knots 1 and 3. In optical cavity it can be semi-glassy mirrors. Signals C 
and D have complex slow amplitude C and D and frequency u. 

Let us assume that near the knot 1 signal in the ring have slow complex amplitude Z n and 
frequency u>. By travelimg (the direction is indicated by arrow) a signal spreads through the 
segment of nonlinear medium of length R\, in which the signal transforms. For example the 
component on doubled frequency, can be excited. Interaction of two these components (in 
the case of its synchonism and by approximation of statiobarity of waves) can be described by 
differential equations 

t— = aa 0, 1— — = pa (2) 
ox ox 

where a(x) is a complex slow amplitude of the component on frequency u, b(x) is a complex 
slow amplitude of a component on frequency 2u, x is a coordinate along a segment of nonlinear 
medium. If on an input to a medium the amplitude of component with frequency u is a(0) ~ Z n , 
then on an exit in some approximation the amplitude of the component of frequency 2u is 
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b(Ri) ~ zl 

Further a signal runs through a filter which passes only component on frequency 2u. After 
that the signal is reflected from a knot 2 and a knot 3, where external pumping D of frequency 
uj is added. On a segment between the third and fourth knot the nonlinear medium of length 
R 2 is located. Transiting through nonlinear medium signal (with a component of frequency 2u 
come from the filter and a component of frequency u added outside) transforms according to 
the equations (2). On an exit from a nonlinear medium the component of a signal on frequency 
uj in some approximation have intensity a{R>i) D*Z^. 

Further the signal transits through the filter which passes component on frequency u. Fi- 
nally, on a knot 1 there is an additional external pumping C on the same frequency. Complex 
amplitudes of external signal and signal in the ring accumulates a(R 2 ) + C. Thus, the signal 
makes the full circle. It have the same frequency u as in the beginning of the path and its 
complex amplitude can be asymptotically wrote as Z n+ \ = F(Z n ) = C + D*Z^. Last expression 
represents complex map. 

Numerical experiment of evolution of a signal in this system is carried out. The equations 
(2) are integrated (by Runge-Kutta method). As a matter of fact the numerical simulation is 
divided into some stage: 1) an integration of a system (2) with parameters a%, 0i and with 
initial conditions a(0) = Z n , 6(0) = gives us a(Ri) and b(Ri); 2) an integration of a system 
(2) with parameters a 2 , 2 and with initial conditions a(0) = D, 6(0) = b(Ri) gives us a(R 2 ) 
and b(R 2 ); 3) accumulation a(R 2 ) + C gives us a new step value of a complex discrete variable 
Z n+ \. At Fig. 3 two charts of a plane of parameter C for an investigated system at different 
values of parameters, characterising the nonlinear media are represented. 

Obviously the Poicare map Z n+ \ = F(Z n ) is not complex analytical (it is not satisfies 
Cauchy-Riemann equations) because system (2) contains non-analytical term, proportional 
to a*. Phenomena of complex analytic dynamics can realise in described model only as a 
approximation. One can see at Fig. 3 "almost ideal" Mandelbrot set on top panel which is 
distorted on bottom panel by the growth of parameters an and f3 2 (leading to amplification of 
" non-analyticity" ) . 

One tool for an analysis of the degree of violation of the complex analyticity, suggested 
in [12] , is the computation of the spectrum of Lyapunov exponents. In particular, for a 
two-dimensional map X n+ i = U(X n ,Y n ), Y n+ i = V(X n ,Y n ) they may be determined via the 
eigenvalues of the matrix 

b = a(X , F )a + (A , Y ) a (X u Y 1 )aL+(X 1 , Y 1 )...a.(X M _ 1 , y M -i)a + (X M _ 1) Y M ^), (3) 

where 

/ dU(X,Y)/dX dU(X,Y)/dY\ 

\dV{X,Y)/dX dV{X,Y)/dY J' U 
In the case of a two-dimensional real map equivalent to an analytic map of one complex variable, 
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Figure 3: Mandelbrot-like set (a, d), which arises for the ring cavity represented at Fig. 2, on 
the chart of the complex plane of a parameter C, responsible to the complex slow amplitude 
of pumping signal C, and Julia-like set (b-c, e-f), which arises at he complex plane of the 
amplitude of the signal Z in the ring in the Poincare section near knot 1. Values of parameters 
are: complex slow amplitude of the pumping signal D D = 3; length of nonlinear elements 
Ri = i?2 = 1-0; parameters, characterising media in nonlinear elements a.\ = 0.01, /3i = 1.0, 
a 2 = 1.0, 2 = 0.01 (arc) and a ± = 0.05, ft = 1.0, a 2 = 1.0, /3 2 = 0.1 (d-f) 
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two Lyapunov exponents must be equal. It may be shown from the Cauchy-Riemann conditions 



dU(X, Y)/dX = dV(X, Y)/dY, 

dU(X, Y)/dY = -dV(X, Y)/dX, 1 ' 

that two eigenvalues coincide at any values of parameters and variables. The same is true for 
the Lyapunov exponents expressed as Ai j2 ~ logAi^. 

The studying ring system possesses four Lyapunov exponents. To compute the Lyapunov 
exponents we used the Benettin algorithm [13]. The procedure consists in simultaneous numer- 
ical solution of the equations (2) and a collection of four exemplars of the linearized equations 
for small perturbations: 

i^- = aa*b + aa*b, i^- = 2/3aa (6) 
ox ox 

with parameters «i, j5\ during x G (0,i?i). and with parameters «2, 02 during x G (0,i?2)- 

Signal transformation in the filters are taken into account as the shift of the variables and the 

perturbation vectors. Pumping signal C gives only impact to the variable. 

At each Poincare section after passing knot 1 in the ring we perform Gram-Schmidt orthogo- 

nalization and normalization for a set of four vectors x J = {Rea , Ima , Re?/, Imfr 7 }, j = 1, 4. 

The Lyapunov exponents are estimated as mean rates of growth or decrease of logarithms of 

the norms of these four vectors: 

1 M 

i=l 

where the norms are evaluated after the orthogonalization but before the normalization. 

Computations show that, depending on the regime, two larger exponents may be negative 
(periodic attractive orbits), positive (chaotic motions) and zero (a border of chaos and quasiperi- 
odic regimes). The second Lyapunov exponent is more or less close to the first one. The other 
two exponents are always negative in the whole domain of existence of bounded dynamical (i.e. 
on the Mandelbrot set) and tend — oo. In the left column of Fig. 4 we present charts of the 
largest Lyapunov exponent on the plane C for the same parameters values as at Fig. 3. Gray 
tones from light to dark correspond to variation of the Lyapunov exponent from to — oo. 
Observe that at central parts of the Mandelbrot set leaves the largest Lyapunov exponent be- 
comes large negative, which corresponds to periodic motions of high stability. At edges of the 
leaves, a thin strip of appearance of positive Lyapunov exponent can happen (chaos). The 
picture is similar to that for the quadratic complex analytic map; see e.g. Ref. [H]. At panel 
(c), distortion of the configuration develops. For example, one can see thick stripe of chaos. 
The leaves lose their round form and separate each other. In the right column of Fig. 4, we 
depict respective charts for the difference of the two larger Lyapunov exponents. The regions 
of large difference of the exponent are visualized by dark gray and black color. It reveals the 
essential deviation from complex analytic dynamics. 
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Figure 4: Charts of the largest Lyapunov exponent (a, c) and for the difference of the two 
larger exponents (b, d) for the ring system at the values of parameters, responsible for nonlinear 
mediae properties: a,\ = 0.01, (3\ = 1.0, = 1.0, 02 = 0.01 (a, b) and a± = 0.05, (3\ = 1.0, 
a 2 = 1.0, (3 2 = 0.1 (c, d). Values of other parameters are: D = 3; Ri = R 2 = 1.0. Uniform 
color means area of unstable dynamics (divergence to infinity). Black color on the diagrams 
corresponds to chaotic dynamics in a bounded domain. More dark tone of gray color in the 
left column corresponds to smaller value of the senior Lyapunov exponent and more strong 
stability of the periodic regime. More dark tone of gray color in the right column means 
larger difference brtween Lyapuniv exponents and therefor more considerable distortion of the 
complex analyticity. 
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